Anomalously large critical regions in power-law random matrix ensembles 
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We investigate numerically the power-law random matrix ensembles. Wavefunctions are fractal up 
to a characteristic length whose logarithm diverges asymmetrically with different exponents, 1 in 
the localized phase and 0.5 in the extended phase. The characteristic length is so anomalously large 
that for macroscopic samples there exists a finite critical region, in which this length is larger than 
the system size. The Green's functions decrease with distance as a power law with an exponent 
related to the correlation dimension. 



PACS number(s): 72.15.Rn, 05.40.-a, 05.45.Df 

The power-law random banded matrix (PRBM) model 
was introduced by Mirlin et al. JjJ to describe a ID sam- 
ple with random long-range hopping. The model at- 
tracted a lot of attention because it can approximately 
represent a variety of different physical systems: an inte- 
grable billiard with a Coulomb scattering center p|, two 
interacting particles in a ID random potential quan- 
tum chaos in a billiard with a non-analytic boundary j| , 
and the Luttinger liquid at finite temperatures Js],^). As 
the model describes a whole family of critical theories, it 
is also interesting for the study of intermediate spectral 
statistics (?) ||]. 

The model is represented by N x N real symmetric ma- 
trices whose entries are randomly drawn from a normal 
distribution with zero mean and a variance depending on 
the distance of the matrix element from the diagonal 
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From field theoretical considerations [gjjg jl0| , |rjj , the 
PRBM model was shown to undergo a sharp transition 
at a = 1 from localized states for a > 1 to delocalized 
states for a < 1. This transition is supposed to be simi- 
lar to an Anderson metal-insulator transition, presenting 
multifractality of eigenfunctions and non-trivial spectral 
compressibility at criticality. The parameter 6 plays a 
role analogous to dimensionality establishing the charac- 
ter of the transition. 

From a perturbative treatment of the non-linear er- 
model and from renormalization group calculations Mir- 
lin considered different regimes depending on the expo- 
nent a. These models are only justified for 6 3> 1, al- 
though the other limiting case can also be studied fol- 
lowing Levitov |l2|. For a > 3/2 the eigenstates are 
localized, but in contrast to usual tight-binding models 
they do not decay exponentially, but with a power-law 
with the same exponent a as the hopping matrix ele- 
ments For 1 < a < 3/2 the states are still local- 
ized, but with a different typical length. In both cases, 
the wavefunctions are expected to have integrable power- 
law tails and the dimensionality to be equal to zero. For 
very large system sizes, the nearest neighbor normalized 



level spacings s should be distributed according to the 
Poisson law, Pp(s) = exp(— s). 

For a < 1, the previous theories predict that all states 
should be delocalized, independently of the value of b. 
The inverse participation ratio (IPR) should then be pro- 
portional to the system size, as for Gaussian orthogonal 
ensembles (GOE), although higher order cumulants of 
the IPR should not necessarily behave in the same way 
as GOE |lo| . For extended states, the nearest neighbor 
level spacing distribution approximately follows Wigner- 
Dyson law 0] , P W d(s) = (tts/2) cxp(-7rs 2 /4) . 

At the critical value a = 1, we expect to have multi- 
fractal eigenstates and critical level statistics, intermedi- 
ate between Poisson and Wigner-Dyson. Recent numer- 
ical calculations at criticality and for 6=1 have shown 
that the nearest level spacing distribution differs from 
the typical one at a metal- insulator transition |H| . 

Our aim in this paper is to explore in detail the char- 
acteristics of the PRBM model in the interesting regime 
around a s» 1 and for the case 6=1 where theoreti- 
cal treatments are not strictly applicable. We do so by 
studying systematically the nearest level spacing distri- 
bution function P(s), the behavior of the Green function 
(GF) as a function of distance and the fractal proper- 
ties of the states. We consider the PRBM Hamiltonian 
given by Eq. (|l]) with the exponent a ranging between 
0.3 and 2. In order to reduce edge effects, we use peri- 
odic boundary conditions, so that our systems are rings 
of length L. We have checked that the results do not 
differ significantly from those obtained for linear chains 
away from the ends. 

We first calculate the eigenvalues of the Hamilto- 
nian matrix by direct numerical diagonalization in or- 
der to study P(s). We consider a small energy window 
(—0.4, 0.4) at the center of the band. For this calcula- 
tion, the system size varies from L = 200 to 5000 and 
the number of random realizations is equal to 10 6 /L. 
We quantitatively analyze the scaling behavior of P(s) 
through the following function of its variance |lq,[l7| 
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which describes the relative deviation of var(s) from the 
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Wigner-Dyson limit due to the finiteness of the system. 
In Eq. (g) var(s) = (s 2 ) - (s) 2 , and 0.286 and 1 are the 
variances of Wigner-Dyson and Poisson distributions, re- 
spectively. We check that the use of other scaling vari- 
ables [ pTs^pToj , instead of 77, do not alter the results. 
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FIG. 1. Scaling variable 77 as a function of a for different 
system sizes: L = 200 (o), 500 (A), 1000 (v) and 3000 (□) 
and 5000 (o). Up to a = 1, 77 is basically independent of L. 

Figure 1 shows the a dependence of 77 for different sys- 
tem sizes: L = 200 (o), 500 (A), 1000 (y) and 3000 (□) 
and 5000 (o). The value of 77 at a = 1 is independent 
of system size and very close to the GOE value, instead 
of being intermediate between Wigner-Dyson and Pois- 
son like in a standard metal-insulator transition. This is 
compatible with the numerical calculations of Ref . [jl5| , 
where it was found that the spectral statistics at a = 1 
differs from that typical of such a transition. For a > 1 
there is a strong increase of 77 with both system size and 
a towards the Poisson limit. For a < 1, 77 is basically 
independent of system size and very similar to the GOE 
value. The small value of rj at criticality makes it diffi- 
cult to decide whether there is a crossover at a = 1, as 
in standard metal-insulator transitions. It also indicates 
that, for 6=1, the present transition is characterized 
by a very large conductance. This characteristic conduc- 
tance would decrease with decreasing b Q . 

Next, we calculate the real space GF G(i,j; E) as the 
inverse of the Hamiltonian matrix. For each value of 
a and L, we obtain G(i,j;0) at the center of the band 
E = for different disorder realizations. We first average 
with respect to the initial position and then we make 
an ensemble average of the logarithm of this quantity, 
obtaining 

F(n) = (ln^|G(M + n;0)|) -InL. (3) 

i 

In this case, the system size varies from L = 300 to 4000 



and the number of random realizations is equal to 10 6 /L. 

In Fig. 2 we plot F(n) — F(Q) as a function of n for 
L = 4000 and for different values of the exponent a = 0.7, 
0.8, 0.85, 0.9, 0.95 and 1. We subtract F(0) in order to 
remove effects associated with variations in the average 
density of states. The central region (tail) of the results 
can be fitted very well by a sum of two power-law terms 
with the same exponent (3 

W^(0) = ln{^ + lr ^}, (4) 

where A is a constant proportional to the density of 
states. Each contribution in Eq. (||) corresponds to prop- 
agation along one arm of the ring. The exponent is 
closely related to the matrix exponent a and weakly de- 
pendent on system size. The fit of Eq. (|J) to the experi- 
mental points is so good that it is indistinguishable from 
the data on the scale used in Fig. 2. Up to the large 
scales used states are not truly extended. The results for 
a > 1 (not shown) are similar to those represented in the 
figure, for a < 1. 
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FIG. 2. F(n) - F(0) as a function of n for L = 4000 and 
a = 0.7, 0.8, 0.85, 0.9, 0.95 and 1, from top to bottom. 

We have tried different extrapolation schemes in order 
to obtain the large scale behavior of the exponent (3. We 
found that a logarithmic dependence on system size L 
fits fairly well the results, as can be seen in Fig. 3 where 
we plot (3 (solid symbols) versus L for a = 0.8, 0.9, 0.95, 
1, 1.05 and 1.1 on a semilogarithmic scale. The data can 
be fitted within errors by straight lines. The value of 
(3 for the critical case, a — 1, is size independent and 
equal to 0.25. For a < 1, j3 decreases with size towards 
0, while for a > 1 it increases towards a. In the limit of 
macroscopic sizes, let us say, L rj 10 s , we have a finite 
regime where the value of the exponent (3 is different from 
zero and so states are not conducting. For infinitely large 
sizes (3 is zero below the transition, 0.25 at the transition 
and a above the transition. 



2 



To clarify the nature of the wavefunctions which give 
rise to this anomalous behavior of the GF, we calculate 
their fractal dimensions in the region around a = 1. We 
obtain the eigenfunctions of the Hamiltonian matrix by 
direct numerical diagonalization for system sizes ranging 
between L = 500 and 10000. For 0.75 < a < 1.25, we 
found that the wavefunctions show a complex structure 
with many sharp peaks of different heights and irregu- 
larly spaced which suggests the possibility of multifrac- 
tal characteristics. We calculate the generalized fractal 
dimensions D q by the standard box counting procedure 
through the expression 



D„ 
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q - 1 (5^0 me) 
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where 5 = 1/ L and is the q-th moment of the prob- 
ability density of the wavefunction in the boxes of size of 
I |2(J. For the evaluation of D\ from the previous equa- 
tion, one has to expand Xq(l) around 5 = 1. In practice 
we compute Eq. (||) at S = 0.1. 
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FIG. 3. Exponent fj characterizing the decrease with dis- 
tance of the GF (solid symbols) and one minus the fractal 
dimension 1 — Di (empty symbols) as a function of L on a 
logarithmic scale for a — 0.8, 0.9, 0.95, 1, 1.05 and 1.1, from 
bottom to top. Solid lines are linear fits to the data. 

We concentrate on the correlation dimension D 2 , which 
reflects the scaling of the density-density correlation func- 
tion. We found that, as for the exponent f3, a logarithmic 
dependence of D 2 with L fits fairly well the numerical 
data. We also found the following empirical relation be- 
tween (3 and D 2 



1 for a < 1 
a for a > 1 
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To show graphically this relationship, we plot in Fig. 
3, together with (3 (solid symbols), the value of 1 — D% 
(empty symbols) as a function of L for a = 0.8, 0.9, 0.95, 



1, 1.05 and 1.1, from bottom to top. In agreement with 
Eq. (^), the data for j3 and 1 — D 2 overlap for a < 1, 
while they are approximately shifted by an amount a — 1 
for a > 1. This shift is consistent with the fact that in 
this regime D 2 tends towards zero, while (3 goes to a. 

At a = 1 , D 2 is practically size independent and equal 
to 0.8. It very slowly tends to for a > 1 and to 1 for a < 
1, but even for (finite) macroscopic sizes remains between 
and 1 for a finite width interval of the parameter a. 
The value of D 2 /d at criticality is larger than for the 
standard Anderson metal- insulator transition in 3D [^0| , 
1.68/3 = 0.56, which indicates that the characteristic 
critical adimensional conductance for the PRBM model 
with b = 1 is larger than for the Anderson model. This is 
in agreement with the small value found for the variance 
of the nearest level spacing at criticality, Fig. 1. 

As for D 2 , the other fractal dimensions are size inde- 
pendent at criticality. We obtained that the information 
(or entropy) dimension D\ is equal to 0.88. Other fractal 
dimensions are £>_ 2 = 1.27, £)_i = 1.14, D 3 = 0.71 
and D4 = 0.65. Trivially the similarity dimension 
Dq = d. In the range of dimensions studied, the de- 
viations with respect to the weak multifractality limit, 
where (D q /d) — 1 oc q, are of the order of 10 %. 

Away from the critical point, the wavefunctions are 
fractal up to a given length, which depends drastically 
on a and diverges as a goes to unity. This length £f rac 
is defined as the system size for which the correlation 
dimension becomes either or 1 in a plot like Fig. 3. In 
Fig. 4 we show log£f rac as a function of a. We see that 
log£f rac diverges at a — 1 in an asymmetric way, unlike 
in standard metal- insulator transitions. In the regime 
a > 1, a perturbative treatment of the non-linear a- 
model predicts the existence of a length scale £ 



log£ 



2a - 1 
'2a- 2 
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which plays the role of the localization length for the 
PRBM model 0,0. Note that £ is not a length charac- 
terizing an exponential decay of the wavefunction, since 
the GF, for system sizes larger than £f rac , decreases with 
distance as a power law with an exponent equal to a. For 
this reason it is the logarithm of £f rac rather than £f rac 
itself what diverges in Eq. (0) with a critical exponent 
and so plays the role of a localization length in a stan- 
dard metal-insulator transition. Eq. (Q) was obtained for 
the regime 6 > 1 and the constant c was estimated to 
be of the order of log b, although one can expect it to be 
qualitatively valid for any value of b. Our data in the 
localized regime can be fitted very well by a curve of the 
form given by Eq. (Q), as can be seen in Fig. 4, where the 
right continuous curve corresponds to this equation with 
c = 1.67. This value of c is larger than predicted by the 
theory if extended to our case 6=1. 

For a < 1, the data of Fig. 4 cannot be fitted by Eq. 
(m), which corresponds to a critical exponent equal to 
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unity, and there is no other theoretical prediction avail- 
able. In this regime the data can be fitted to a curve of 
the form 



log£ 



a-1 
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with the critical exponent v equal to 0.5. The left solid 
line in Fig. 4 corresponds to Eq. (^J) with v = 0.5 and 
the fitting parameter c = 1.92. It is the first time to our 
knowledge that a metal-insulator transition is approach 
by two different critical exponents, which is in contradic- 
tion with a single parameter scaling theory. The possi- 
bility of the existence of two characteristic lengths has 
been already pointed out in a different one-dimensional 
model Q . We note that we are able to obtain two expo- 
nents numerically because the transition point is known 
exactly (a — 1), which constitutes a key advantage of 
this model. 
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FIG. 4. Length £f rac as a function of a. Solid lines corre- 
spond to Eq. (0) with c = 1.67, and Eq. (|) with c = 1.92 and 
v = 0.5. 

To summarize, we have explored the critical region of 
the PRBM model in the case 6=1 inaccessible to field 
theoretical methods. Our results, although compatible 
with theoretical predictions, show new light on the spe- 
cific form of the GF and on the multifractal structure 
of the wavefunctions near the critical point. The wave- 
functions are fractal up to a given length which diverges 
asymmetrically with two different critical exponents at 
a = 1. It is the first time that such an asymmetry has 
been obtained in part due to the exact knowledge of the 
position of the critical point. The possibility of a similar 
asymmetric behavior in other metal-insulator transitions 
should be reconsidered. We note that a fitting procedure 
of the raw data without an exact knowledge of the criti- 
cal value of the transition parameter is compatible with 
equal exponents in the diverging lengths, above and be- 
low the transition. The fractal behavior is correlated with 



a power law decay of the GF with small exponents. The 
sum of the exponent characterizing the GF decay and the 
correlation dimension is approximately equal to unity in 
the localized phase and to a in the extended phase. All 
this should have important consequences in the problem 
of two interacting particles and in the other problems re- 
lated to the PRBM model although their study is clearly 
outside the scope of this work. 
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